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ABSTRACT 


The  beam  element  available  in  computer  program  ADINA  is 
inadequate  for  transient  analysis  of  eccentrically  stiffened 
shell  structures,  particularly  when  the  lateral  stability  of  the 
stiffener  is  of  concern.  As  an  alternative  to  modeling  shell 
stiffeners  with  numerous  continuum  or  transition  elements,  a 
stiffener  modeling  technique  based  on  the  multipoint  constraint 
option  in  ADINA  is  presented.  This  technique  leads  to  sig¬ 
nificant  reductions  in  the  number  of  elements  and  solution 
degrees  of  freedom  needed  for  accurate  stiffener  modeling,  yet 
allows  inclusion  in  the  analysis  of  effects  of  out-of-plane 
web  distortion,  longitudinal  warping,  and  torsion.  Stiffeners 
having  various  cross  sectional  geometries  and  boundary  conditions 
have  been  modeled,  and  predicted  response  correlates  well  with 
experimental  data.  The  approach  is  of  practical  significance 
for  large  stiffened  shell  problems,  especially  nonlinear  analysis. 

ADMINISTRATIVE  INFORMATION 

The  work  described  in  this  report  was  funded  by  the  Naval  Sea  Systems  Command 
(SEA  3221)  and  was  performed  under  David  W.  Taylor  Naval  Ship  Research  and  Develop¬ 
ment  Center  Program  Element  Submarine  Combat  Protection,  Task  Area  SF  434  00  001 
and  Work  Unit  1750-294. 

INTRODUCTION 

Transient  analysis  of  eccentrically  stiffened,  shell  structures  can  be  quite 
complex.  Under  severe  loading  conditions  such  as  impacts  or  blast  waves,  shell 
stiffeners  can  become  laterally  unstable.  Lateral  instability  of  a  stiffener  leads 
to  a  sudden  loss  of  strength  of  transverse  sections  and  can  result  in  catastrophic 
structural  failure.  The  significance  of  lateral  instability  of  stiffeners  as  a 
primary  mode  of  structural  failure  is  indicated  in  Figure  1. 

Combinations  of  inplane  and  lateral  loading  can  produce  significant  out-of¬ 
plane  deformations  in  a  shell  stiffener.  Unfortunately,  lateral  instability  of 

stiffeners  caused  by  dynamic  loading  is  poorly  understood.  As  pointed  out  by 
1* 

Adamchak,  the  primary  difficulty  encountered  in  attempts  to  analyze  lateral  be¬ 
havior  of  shell  stiffeners  is  the  inability  to  properly  account  for  web  deformations 
or  to  adequately  model  complex,  higher  order,  out-of-plane  modal  displacement 
patterns.  To  include  the  effects  of  these  phenomena  in  a  general,  nonprismatic 
"beam"  element  would  be  a  formidable  task  and  is  probably  beyond  the  scope  of  the 
state  of  the  art  of  finite  element  technology.  As  an  example,  the  three-dimensional 

*A  complete  listing  of  references  is  given  on  page  21. 
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beam  element  in  ADINA  (a  finite  element  program  for  automatic  dynamic  incremental 
2  3 

nonlinear  analysis)  9  assumes  a  prismatic  cross  section  with  no  centroidal  offsets. 

Although  the  beam  element  of  ADINA  does  provide  torsional  stiffness  about  the 

principal  axis,  it  cannot  account  for  out-of -plane  deformation  and  warping  of  the 

cross  section.  An  improved  isoparametric  curved  beam  element  has  been  formulated  by 

4 

Bathe  and  Bolourchi.  This  element  is  fully  compatible  with  both  the  three- 
dimensional  and  thin-shell  elements.*  Unfortunately,  it  lacks  the  needed  generality 
for  modeling  arbitrary  stiffener  geometries  and  out-of-plane  deformation  of  the 
cross  section. 

An  alternate  modeling  strategy  that  can  overcome  previously  enumerated  dif¬ 
ficulties  is  to  fully  model  the  stiffener  as  an  assemblage  of  continuum  or  thin- 
shell  elements.  Unfortunately  for  arbitrary  stiffener  geometries,  such  an  approach 
can  require  numerous  elements  and  excessive  degrees  of  freedom  for  accurate  modeling. 

In  this  report,  a  constrained-displacement  modeling  technique  is  presented. 

The  ADINA  displacement  constraint  option,  together  with  three-dimensional  solid  or 
thin-shell  elements,  is  used  to  model  shell  stiffeners  of  arbitrary  geometry.  The 
number  of  elements  and  solution  degrees  of  freedom  can  be  substantially  reduced,  even 
though  effects  of  out-of-plane  deformations  are  retained  in  the  analysis. 

STIFFENER  DEFORMATION  UNDER  COMBINED  LOADINGS 
The  deformed  geometry  of  a  flexible  shell  stiffener  under  the  action  of  in¬ 
plane  and  lateral  loading  is  shown  in  Figure  2.  Lateral  deformations  can  be 
characterized  as  a  combination  of 

1.  Out-of-plane  bending  and  distortion  of  the  stiffener  (u,v) 

2.  Longitudinal  warping  or  twisting  of  the  cross  section  (angle  S) 

3.  Rotation  and  bending  of  the  shell  plating  to  which  the  stiffener  is 
attached  (angle  a). 

This  report  is  mainly  concerned  with  curved  shell  structures  and  combinations  of 
asymmetric  distributed  pressures  p(x)  and  eccentric  concentrated  forces  f^  (Figure 
3)  which  load  the  stiffener  primarily  in  transverse  compression  and  torsion.  These 
loads  pose  the  most  significant  threat  to  lateral  stability  of  the  stiffener  because 
they  tend  to  lessen  transverse  and  torsional  resistance  of  the  stiffened  shell 


*It  is  ur  ^<*tood  .at  the  isoparametric,  curved  beam  element  will  be  available 
in  the  1981  rele**'  oi  aDINA. 
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section  and,  also,  cause  the  stiffener  to  distort  away  from  the  stable  equilibrium 
configuration.  Consequently,  an  adequate  model  of  a  stiffened  shell  with  arbitrary 
cross  section  geometry  should  include  sufficient  detail — node  and  gaussian- 
integration  points — to  map  complex  lateral  distortion  patterns  as  well  as  to  account 
for  large  gradients  in  the  stress  and  strain  fields. 

CONSTRAINED-DISPLACEMENT  STIFFENER  MODEL 

Some  examples  of  finite  element  modeling  for  flexible  shell  stiffeners  are  shown 
in  Figure  4.  The  continuum  element  model  is  fully  compatible.  In  contrast,  the 
ADINA  isoparametric,  thin-shell  element  formulation  allows  only  five  degrees  of 
freedom — three  translations  and  two  rotations — at  a  node  and  requires  the  midsurface 
normal  vector  to  be  specified  at  each  midsurface  node.  This  requirement  leads  to 
a  discontinuity  at  shell  intersections  which  must  be  resolved  to  prevent  a  singular 
stiffness  matrix.  In  Reference  2,  it  is  suggested  that  such  discontinuities  can  be 
overcome  with  constraint  equations  of  the  type  indicated  in  Figure  5.  It  is  proposed 
herein  that  the  ADINA  linear,  multipoint  displacement  constraint  option  be  used  in  a 
more  general  fashion  to  overcome  the  previously  outlined  difficulties  and,  also,  to 
reduce  the  number  of  elements  and  solution  degrees  of  freedom  required  for  more 
effective  analysis  when  a  continuum  element  model  is  preferred. 

In  ADINA,  linear,  multipoint  displacement  constraint  equations  have  the  general 

form 

N 

u  .  =  y  a,  u  (I) 

ij  X  v  k  mn 

k=l 


where  k  =  1 ,  2 ,  3 . . . N 

1  £  m  NUMNP  (number  of  node  points) 
n  -  1,  2,  or  3 


In  Equation  (1),  n  is  the  number  of  independent  translational  or  rotational 

degrees  of  freedom,  a  are  weighting  factors,  and  u  is  the  displacement  or  ro- 

k.  ij 

tational  degree  of  freedom  at  the  i  th  node  point  in  or  about  the  j  th  direction  or 


axis . 
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If  a  master-slave  nodal  dependency  is  established  for  common  nodes  at  inter¬ 
sections  of  shell  and  stiffener  or  web  and  flange,  then  stiffened  shell  sections 
(Figure  6)  can  be  easily  modeled.  Development  of  suitable  constraint  equations  for 
the  dependent  degrees  of  freedom  is  straightforward.  As  an  example,  consider  the 
continuum  element  model  in  Figure  7.  To  a  first  order  approximation,  the  angle  of 
rotation  of  the  intersection  of  shell  and  stiffener  about  the  normal  to  the  cross 
sectional  plane  is  defined  as 


a 

z 


(2) 


and  u,  v,  and  w  are  components  of  displacement  in  the  x,  y,  and  z  or  r,  s,  and  t 
coordinate  directions.  If  it  is  assumed  that  the  angle  of  intersection  remains 
constant  during  deformation,  then  the  dependent  degrees  of  freedom  become 


w 

u-  =  +  —  a  =  u0  -  a  v0  +  a  v, 

1  3  2  z  3  o3  o  4 


(3) 


V1  v3 


(4) 


W1  w3 
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w 
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w  *  w 

2  3 
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(8) 


where  a 


•*ft) 


(9) 


Equations  (3) -(9)  are  valid  in  either  global  Cartesian  or  local  skewed  co¬ 
ordinate  reference  systems  as  long  as  the  normal  to  the  cross  sectional  plane  is  in 
the  direction  of  a  principal  axis  because  displacements  of  interest  will  always  lie 
in  a  principal  plane.  If,  however,  the  normal  does  not  coincide  with  a  principal 
direction,  for  example,  for  a  curved  shell  for  which  displacements  are  calculated 
in  terms  of  global  Cartesian  coordinates  (Figure  8),  the  components  of  the  rotate;  ;1 
displacement  in  the  principal  directions  are  obtained  with  a  simple  transformation. 
For  such  cases,  the  angle  of  rotation  about  the  normal  is  defined  as 


V4  “  v3 

t  =  - 

t  t 


and  the  increment  of  displacement  due  to  this  rotation  becomes 


Ar "  f  \  "  i  (r)  W 


The  components  of  displacement  in  the  principal  directions  are 


Ax  =  sin  0  Ar 


Az  =  cos  0  Ar 


and  constraint  equations  for  the  dependent  degrees  of  freedom  become 


=  u3  +  Ax  «  u3  -  a^v^  +  a^v4 


V1  '  v3 


+  Az  «  +  a^v^  -  a^v^ 


1 


i 


u«  *  u~  -  Ax  *  uu  +  a,v~  -  a, v , 
2  3  3  1  3  1  4 


V2  “  V3 


(17) 

(18) 


w2  =  w3  -  Ae  =  w3  -  a2v3  +  a2v4 


(19) 


sin  0  (20) 


cos  0  (21) 

EXAMPLE:  LATERAL  STABILITY  OF  A  STIFFENER  UNDER  IMPACT  LOADING 
The  constrained-displacement  modeling  technique  was  used  to  investigate  the 
stiffener  lateral  stability  of  a  stiffened  cylindrical  shell  subjected  to  a  known 
impact  loading.  Figure  9  shows  the  shell  geometry  and  loading  conditions.  Figure  10 
shows  the  corresponding  ADINA  finite  element  model.  Three-dimensional  continuum 
elements  were  used  to  model  both  the  shell  plating  and  the  stiffeners.  Only  one 
T-stiffener  was  fully  modeled.  The  other  was  modeled  as  a  dynamically  equivalent 
bar  stiffener.  A  90-degree  sector  of  the  shell  was  modeled  with  symmetry  imposed  on 
nodal  displacements  along  the  longitudinal  boundaries.  Transverse  boundaries  were 
modeled  as  pinned  supports.  Displacement  constraints  were  imposed  at  intersections 
of  shell  and  stiffener  and  web  and  flange. 

In  Figure  11,  predicted  radial  response  of  the  T-stiffener  is  compared  to 
experimentally  measured  response  data.  Correlation  between  the  two  is  quite  good. 
Lateral  response  of  the  flange  is  shown  in  Figure  12.  The  higher  order  modal  dis¬ 
placement  pattern,  characteristic  of  out-of-plane  stiffener  deformation,  is  clearly 
indicated.  In  the  study,  it  was  observed  that  the  peak  lateral  response  of  the 
frame  was  30  percent  of  the  peak  radial  response,  and  peak  stresses  in  the  flange 
were  comparable  to  peak  stresses  in  the  hull  plating  in  the  vicinity  of  the  impact 
point.  In  contrast,  the  bar  stiffener  model  imparted  a  higher  degree  of  torsional 
resistance  to  the  shell  plating  and  grossly  underpredicted  the  stiffener  stresses. 
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These  results  demonstrate  the  accuracy  of  the  constrained-displacement  modeling 
technique  and  provide  insight  into  the  practical  significance  of  out-of-plane 
stiffener  deformations. 

CONCLUSIONS 

A  constrained-displacement  modeling  technique  for  eccentrically  stiffened  shell 
structures  has  been  presented  which  allows  the  effects  of  out-of-plane  deformations 
of  the  stiffener  to  be  included  in  an  analysis.  For  the  transient  analysis  of  an 
impulsively  loaded  shell  structure,  predicted  response  correlates  well  with  experi¬ 
mentally  measured  data.  This  technique  overcomes  the  difficulties  encountered  by 
using  generalized  beam  elements  while  minimizing  the  number  of  continuum  elements 
and  solution  degrees  of  freedom  required  for  accurate  stiffener  modeling. 


Figure  1  -  Example  of  Stiffener  Lateral  Instability 
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Figure  3a  -  Asymmetric  Pressures 


Figure  3b  -  Eccentric  Forces 
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Figure  3c  -  Equivalent  Stiffener  Loads 


Figure  3  -  Significant  Shell  Stiffener  Loading 
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18  NODES 
5  ELEMENTS 
54  DEGREES  OF 
FREEDOM 


Figure  4a  -  Continuum  Model 


Figure  4b  -  Thin  Shell  Model 


Figure  4  -  Stiffener  Modeling  Techniques 


NODES  i  AND  j  HAVE  SAME 
COORDINATES  AND  DE¬ 
GREES  OF  FREEDOM 


Figure  5  -  Modeling  Techniques  for  Shell  Intersections 
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Figure  6a  -  Continuum  Model 


THICK  SHELL 


Figure  6b  -  Hybrid  Model 

Figure  6  -  Constrained-Displacement  Stiffener  Models 
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Deformations  at  a  Shell  Stiffener  Intersection 
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Figure  9a  -  Shell  Geometry 
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Figure  9b  -  Impact  Loading 


Figure  9  -  Stiffened  Cylindrical  Shell  under  Impact  Loading 


Figure  10a  -  Shell  Plating  Mesh 


T  BAR 

STIFFENER  STIFFENER 

Figure  10b  -  Stiffener  Modeling 


Figure  10  -  Model  for  ADINA  Analysis  of  Stiffener  Lateral  Stability 
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Figure  11  -  Radial  Response  of  Stiffener 
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Figure  12a  -  Response  Predicted  by  ADINA 


Figure  12  -  Lateral  Response  of  Stiffener 
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25 

175.2 

T.  Giacofci 

1 

175.2 

S.  Zilliacus 

1 

175.2 

P.N.  Roth 

1 

175.2 

F.  Fortin 

1 

175.2 

S.  Walter 

1 

175.2 

M.  Hoffman 

1 

175.2 

T.  Tor id is 

1 

175.3 

W.  Conley 

1 

177 

1 

177.1 

V.  Bloodgood 

1 

177.1 

F.  Costanzo 

1 

177.7 

(m) 

10 

5211.1 

Reports  Distr 

1 

522.1 

Unclassified 

1 

522.2 

Unclassified 

